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ABSTRACT 

We present the results of M31 globular cluster (GC) X-ray source survey, based on the data of XMM- 
Newton and Chandra observations covering ~ 6100 arcmin^ of M31. We detected 43 X-ray sources 
coincident with globular cluster candidates from various optical surveys. 

The inferred isotropic X-ray luminosities of GC sources lie between ~ 10'^^ and ~ 10'^^ erg s^^ in the 
0.3 — 10 keV energy band. The spectral properties of 31 brightest sources from our sample were found to 
be similar to that of the low mass X-ray binaries located in the bulge and globular clusters of the Milky 
Way Galaxy. The spectral distribution of M31 GC X-ray sources is consistent with that derived for the 
bulge of M31 and other nearby galaxies of different morphological type. Several sources demonstrate 
a correlation between the level of X-ray flux and hardness of their energy spectrum reminiscent of the 
Galactic Z and atoll sources. 

We found that ^ 80% of the M31 GC sources with multiple flux measurements available show signifi- 
cant variability on a time scales from days to years. The X-ray source RX J0043. 2-1-4127, coincident with 
GC Bo 163, has been found to show recurrent transient outbursts with peak luminosities of ^ lO'^^ ergs 
s~^. Several sources in our sample show significant variability on a time scale of individual observations, 
ranging from aperiodic fluctuations to regular dipping. 

The X-ray luminosity function of GC sources is found to be significantly different from that of the point 
sources in the bulge and disk of M31. The luminosity distribution of M31 GC sources has ~ 10 times 
higher peak luminosity and much higher fraction of bright sources than the Milky Way GC distribution. 
Six persistent sources in our sample (or 14% of the total number) have luminosities exceeding 10^^ 
ergs s~^ during all observations, and three other sources occasionally exceed that luminosity level. Our 
observations indicate that GC sources make dominant contribution to the bright source counts in the 
areas of M31 covered by the survey: ~ 40% of the total number of sources with luminosities above 10^^ 
ergs s~^ reside in GCs with fraction of GC sources rising to 67 — 90% for the luminosities above 10"^* ergs 
s'-^. The contribution of the GC sources to the total number of bright sources found in M31 is much 
higher than in the Milky Way galaxy, but surprisingly close to the early- type galaxies. We found that 
brightest M31 GC sources tend to reside at large galactocentric distances outside the central bulge. 

We found that globular clusters hosting bright X-ray sources are optically brighter and more metal 
rich than the rest of M31 globular clusters, in agreement with previous studies. The brightest sources 
with luminosities above ~ 10'^* ergs s~^ show tendency to reside in more metal poor clusters. 

The remarkable similarities between the properties of the M31 GC X-ray sources and that of the 
Galactic neutron star LMXBs allow us expect most of the persistent M31 GC X-ray sources to be 
LMXB systems with neutron star primaries. However, the current X-ray spectral and timing data can 
not rule out the possibility of finding an active accreting black holes in our GC source sample. 



1. INTRODUCTION 

Globular clusters (GCs) provide us with crucial informa- 
tion on galaxy structure and formation mechanisms. They 
are also ideal laboratories for studying stellar populations 
and evolution of dense stellar systems. X-ray surveys re- 
vealed a number of bright X-ray sources associated with 
Milky Way (MW) GCs and identified with low-mass X- 
ray binaries (LMXBs) (Hertz & Grindlay 1983, Verbunt 
et al. 1995). The ratio of LMXBs to stellar mass is two 
orders of magnitude higher for GCs than for the rest of our 
Galaxy (Liu et al. 2001). This overabundance of LMXB 
has been explained by dynamical effects such as tidal cap- 
ture and three/four body interactions in the high stellar 
density environment (Clark 1975, Fabian, Pringle & Rees 
1975). The number of known Galactic GCs hosting bright 



X-ray sources {Lx > lO'^^ ergs s^^) is pretty small: there 
are fourteen such systems (Liu et al. 2001). The pres- 
ence of Type I X-ray bursts indicates that all luminous 
GC sources contain accreting neutron stars rather than 
black holes. With the advent of the modern X-ray obser- 
vatories, it has become possible to study the X-ray proper- 
ties of GCs associated with the nearby galaxies (Angelini, 
Loewenstein & Mushotzky 2001; Di Stefano et al. 2002; 
Kundu, Maccarone, & Zepf 2002). Based on a large num- 
ber of objects, the extragalactic GC X-ray sources allow 
a broad range of statistical population studies, impossible 
with a limited sample of Galactic GC X-ray sources. 

The Andromeda Galaxy (M31), the closest giant spiral 
galaxy to our own, is a unique object for the study of op- 
tical and X-ray astronomy. Its proximity and favorable 
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orientation allow to observe stellar populations over the 
full extent of the galaxy at a nearly uniform distance, and 
with less severe effects of line-of-sight contamination from 
interstellar gas and dust. Due to similarities between the 
two galaxies, the results from the study of M31 provide an 
important benchmark for comparison with the results from 
the study of our own Milky Way Galaxy. The Andromeda 
Galaxy possesses a significantly more populous globular 
cluster system than the Milky Way with more than 430 
confirmed candidate members (Barmby et al. 2000). Op- 
tical studies of the Milky Way and M31 globular cluster 
systems have revealed that these two populations exhibit 
some remarkable similarities (van den Bergh 2000). 

M31 was observed extensively with Einstein, ROSAT, 
Chandra and XMM missions, detected hundreds of sources, 
with bright GC X-ray sources among them (Trinchieri & 
Fabbiano 1991; Primini et al. 1993; Supper et al. 1997; 
Shirey et al. 2001). Twenty-one X-ray sources discov- 
ered with Einstein were tentatively identified with GCs. 
Based on these early results, it has been noted that the 
luminosities of the M31 GC sources are higher than any 
observed for Galactic GCs (Long & van Speybroek 1983), 
and that the fraction of M31 GCs with X-ray sources might 
be larger than the fraction in the Milky Way (Battistini et 
al. 1987). Observations of the central bulge of M31 with 
ROSAT/HRl detected 18 GC X-ray sources and placed up- 
per limits on the X-ray flux from 32 other globular clusters 
(Primini et al. 1993). Primini, Forman & Jones (1993) 
found that the luminosity distribution of M31 GC sources 
contains a larger population of high luminosity sources 
than that of the Milky Way GC X-ray sources, although 
the peak luminosities of the two distributions were roughly 
comparable. The ROSAT/PSPC survey covering most of 
M31 revealed 31 X-ray sources coincident with M31 GC 
candidates (Supper et al. 1997, 2001). Based on this larger 
sample, Supper et al. (1997, 2001) concluded that there 
is no significant difference between the luminosity distri- 
butions of the M31 GC X-ray sources and their Galactic 
counterparts. Recent Chandra observations covering cen- 
tral bulge of and two regions in the northern and southern 
disk revealed 28 GC X-ray sources with 15 of them newly 
discovered (Di Stefano et al. 2002). Approximately ~ 30% 
of all GC sources detected in Chandra observations have 
0.5 — 7 keV luminosities above 10^^ ergs s~^ and ~ 10% 
of all the sources have luminosities of 10^^ ergs s~^ and 
higher. Both the peak X-ray luminosity and relative frac- 
tion of the bright sources with Lx > 10^^ ergs s~^ were 
found to be significantly higher for the Chandra GC source 
sample than for the GC sources in our Galaxy (Di Stefano 
et al. 2002) in contrast to the earlier i? 05^7 results (Sup- 
per et al. 1997, 2001). The limited sensitivity of Chan- 
dra observations allowed to study spectral properties and 
variability of only five most luminous sources found to be 
consistent with properties of the Galactic sources. 

We present the results of M31 globular cluster (GC) X- 
ray source survey, based on the data of XMM-Newton and 
Chandra observations covering ~ 6100 arcmin^ of the An- 
dromeda galaxy (M31). The unique combination of the 
unprecedented sensitivity of XMM-Newton and superior 
spatial resolution of Chandra observations allows accu- 

■^See http:/ /xmm.vilspa.csa.cs/user 
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rate localization and detailed study of spectral and timing 
properties of the GC X-ray sources in both central regions 
of M31 and at large galactocentric distances. We concen- 
trate on the analysis of spectral properties and variability 
of individual bright sources and study a group properties 
of all GC X-ray sources detected in our survey. We study 
the relation between the properties of the X-ray sources 
and that of the globular clusters hosting them. Finally, 
the properties of M31 GC X-ray sources are compared to 
the properties of low- mass X-ray binaries LMXB in M31, 
Milky Way and other galaxies of different morphological 
types. 

2. OBSERVATIONS AND DATA ANALYSIS 

In the following analysis we used the data of 10 XMM- 
Newton observations of various M31 fields (Table 1; Fig. 
1). All XMM-Newton observations described in this paper 
were performed as a part of the Performance Verification 
and Guaranteed Time Programs (Shirey et al. 2001, Tru- 
dolyubov et al. 2002a, Trudolyubov et al. 2002b). Four 
XMM observations cover the central bulge of M31, five - 
its northern and southern disk regions, and one - a distant 
M31 globular cluster Gl (Mayall II) field. In the following 
analysis we use the data of three European Photon Imag- 
ing Camera (EPIC) instruments: two EPIC MOS detec- 
tors (Turner et al. 2001) and EPIC-pn detector (Strueder 
et al. 2001). During all observations EPIC instruments 
were operated in the full window mode (30' FOV) with 
medium or thin optical blocking filter. 

We reduced EPIC data with the XMM-Newton Science 
Analysis System (SAS v 5.3) ^. We performed standard 
screening of the EPIC data to exclude time intervals with 
high background levels. Images in the celestial coordi- 
nates with a pixel size of 2" have been accumulated in 
the 0.3 - 7.0 keV energy band for the EPIC-MOSl, M0S2 
and pn detectors. To generate lightcurves and spectra of 
X-ray sources, we used elliptical extraction regions with 
semi-axes size of ~ 20 — 80" (depending on the distance 
of the source from the telescope axis) and subtracted as 
background the spectrum of adjacent source-free regions 
with subsequent normalization by a ratio of the detector 
areas. We corrected the count rates of the sources for the 
vignetting of the XMM telescope, based on the Current 
Cahbration Files. We used data in the 0.3 — 10 keV energy 
band because of the uncertainties in the calibration of the 
EPIC instruments outside this range. All fluxes and lumi- 
nosities derived from spectral analysis apply to this band, 
unless specified otherwise. We used spectral response ma- 
trices generated by XMM SAS tasks. The EPIC count 
rates were converted into energy fiuxes in the 0.3 — 10 keV 
energy band using analytical fits to the spectra for brighter 
sources or PIMMS^, assuming standard parameters: an 
absorbed simple power law model with Nu = 7 x 10^° 
cm~^ and photon index a = 1.7 (Shirey et al. 2001). 

We analyzed a series of publicly available Chandra data 
for M31 fields containing optically identified GC candi- 
dates consisting of 36 observations with the Advanced 
CCD Imaging Spectrometer (ACIS) and 3 High Resolution 
Camera (HRC) observations (Table 2; Fig. 1). A more 
detailed description of these observations can be found in 
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Garcia et al. (2000), DiStefano ct al. (2002), Kong ct al. 
(2002) and Kaarct ct al. (2002). 

The data of Chandra observations was processed using 
the CI AO v3.0^ threads. We performed standard screen- 
ing of the Chandra data to exclude time intervals with 
high background levels. Images in the celestial coordi- 
nates with a pixel size of 1" have been accumulated in the 
0.3 — 7.0 keV energy band for source detection procedure. 
Only data in the 0.3 — 7.0 and 0.5 — 7.0 keV energy ranges 
was used in the spectral analysis of ACIS-S and ACIS-I 
observations. To estimate energy fluxes and spectra, we 
extracted counts within source elliptical region with semi- 
axes of 3"— 40" (depending on the source distance from the 
telescope axis). Background counts were extracted from 
the adjacent source-free regions with subsequent normal- 
ization by ratio of detector areas. To account for the con- 
tinuous degradation in the effective low-energy quantum 
efficiency of the ACIS detectors^, we performed correction 
of the spectral response files applying ACISABS model to 
the ARF files. For a few brightest GC sources pilcup frac- 
tion in their ACIS spectra is estimated to be at the level 
of ~ 10 — 30%. In some cases, to correct for a pileup effect 
in the ACIS spectra, we extracted source spectrum from 
the PSF wings, where the source count rate is low enough 
that pilcup effect can be ignored. We also used Sherpa^ 
implementation of the pileup model by Davis (2001,2003). 
The ACIS count rates were converted into energy fluxes in 
the 0.3 — 10 keV energy band using analytical fits to the 
spectra for brighter sources or using PIMMS with spec- 
tral parameters derived from more sensitive observations 
or standard parameters otherwise. For the HRC count 
rates the source fluxes were estimated with PIMMS, based 
on the observed count rates and spectral parameters mea- 
sured with AMM/EPIC and Chandra/ ACIS. 

For both XMM-Newton and Chandra data, X-ray 
sources were detected with the program based on the 
wavelet decomposition algorithm, set at a Aa threshold. 
For our current analysis of XMM-Newton observations, 
we expect error in the source position determination to 
be dominated by residual systematic error of the order 
1 - 3". 

We performed a spectral and timing study of the bright- 
est X-ray sources with total number of counts larger than 
300. Spectra were grouped to contain a minimum of 20 
counts per spectral bin and fit to analytical models using 
the XSPEC v.ll^ fitting package (Arnaud 1996). EPIC- 
pn, MOSl and M0S2 data were fitted simultaneously, but 
with independent normalizations. 

We studied timing properties of all bright X-ray sources 
detected in our observations. Fourier power density spec- 
tra (PDS) were produced using the lightcurves in the 
0.3 — 7.0 keV energy band. Then we performed folding 
analysis in the vicinity of the frequency peaks identified 
from PDS. We used standard XANADU/XRONOS v.5 1° 
tasks to perform analysis of the timing properties of bright 
X-ray sources. 

In the following analysis we assume a source distance of 

®http:/ /asc. harvard.edu/ciao/ 

^http:/ /cxc. harvard.edu/cal/Acis/Cal prods/qcDcg/ 
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760 kpc (van den Bergh 2000). 

3. M31 GC X-RAY SOURCES: IDENTIFICATION. 

The positions of 43 X-ray sources detected in the XMM- 
Newton and Chandra observations of M31 are consistent 
with globular cluster candidates from optical surveys. We 
used the following catalogs for identifying GC X-ray source 
candidates: the Bologna catalog (Battistini et al. 1987), 
the catalog by Magnier (1993), and the HST globular clus- 
ter candidate catalog (Barmby & Huchra 2001) - with 
search radius of 3". The expect that 4 of our GC X-ray 
candidates could be spurious matches. The information 
on the properties of the X-ray sources and their optical 
identifications is shown in Table 3. 

The X-ray images of M31 with GC X-ray source posi- 
tions circled are shown in Fig. 2 and 3. The X-ray image 
shown in Fig. 2 was constructed combining the data of 
the XMM-Newton/EPlC-MOSl, M0S2 and pn cameras 
in the 0.3 — 7.0 keV energy band. The X-ray image shown 
in Fig. 3 combines the data of Chandra/ ACIS and HRC 
detectors. 

4. NOTE ON THE EFFECT OF MULTIPLE UNRESOLVED 
X-RAY SOURCES 

Recent Chandra observation of the Galactic globular 
cluster M15 has shown two distinct luminous X-ray sources 
separated by less than 3" (White & Angelini 2001). The 
spatial resolution of both Chandra (^ 0.5") and XMM 
(~ 10") is not sufficient to resolve most of possible mul- 
tiple X-ray sources within one cluster in M31. Therefore, 
some of the M31 GC X-ray sources from our sample can 
be unresolved composites of separate sources (Angelini, 
Loewenstein & Mushotzky 2001; Di Stefano et al. 2002). 
It should thus be noted, that the spectral blending and 
superposition of variability of individual sources may com- 
plicate direct comparisons with Galactic globular cluster 
sources. 

5. SPECTRAL PROPERTIES OF THE M31 GLOBULAR 
CLUSTER X-RAY SOURCES 

Twenty seven and twenty three sources detected with 
XMM-Newton and Chandra have a sufficient number of 
counts to allow crude or detailed spectral modeling. Com- 
bining the above two samples gives us a high-quality infor- 
mation on the spectral properties of 31 bright GC X-ray 
sources. The representative spectra of bright GC sources 
obtained with XMM-Newton/'EPIC and Chandra/ ACIS 
are shown in Fig. 4. 

The spectra of globular cluster candidates were fitted 
with a variety of spectral models using XSPEC vll. We 
first considered simple one-component spectral model: an 
absorbed simple power law. The results of fitting this 
model to the source spectra are given in Tables 4 and 5. 
The spectra of all these objects can be generally described 
by an absorbed simple power law model with photon in- 
dex of ~ 0.8 — 2.8 and an equivalent absorbing column of 
~ (0.4-7.5) X 10^^ cm-2 (Table 4,5). The corresponding 
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isotropic luminosities of the globular cluster sources differ 
by four orders of magnitude and fall between ^ 10'*''' and 

10^^ ergs s~^ in the 0.3— 10 keV energy band, assuming 
a distance of 760 kpc. 

In most cases, we obtained acceptable fits using a power 
law spectral model (Table 4,5). However, for almost 
all brighter sources, a complex spectral models are re- 
quired. For several sources with high luminosities (sources 
##22,26,32,42,43 in Table 3), the models with quasi- 
exponential cut-off at ^ 2 — 8 keV or two-component mod- 
els describe the energy spectra significantly better than a 
simple power law. To approximate the spectra of these 
sources, we used an absorbed power law model with ex- 
ponential cut-off (XSPEC CUTOFFPL model), a Comp- 
tonization model and a two-component models. 

For the Comptonization model approximation, we used 
the XSPEC model COMPTT (Sunyaev & Titarchuk 
1980, Titarchuk 1994, Titarchuk & Lyubarskij 1995). This 
model includes self-consistent calculation of the spectrum 
produced by the Comptonization of the soft photons in 
a hot plasma. It contains as free parameters the tem- 
perature of the Comptonizing electrons, kT^, the plasma 
optical depth with respect to the electron scattering, t 
and the temperature of the input Wien soft photon dis- 
tribution, kTo. A spherical geometry was assumed for the 
Comptonizing region. 

A cut-off power law and a Comptonization models gave 
reasonably good fit to the spectra of all sources consid- 
ered, and the results are shown in Table 6. The model pa- 
rameters derived for the M31 GC X-ray sources resemble 
those observed for bright Galactic low mass X-ray binaries 
(White, Stella & Parmar 1988; Christian & Swank 1997; 
Church & Balucinska-Church 2001). 

The spectra of many luminous Galactic LMXB are well 
fit with a two-component model consisting of a soft black 
body-like component which might represent emission from 
an optically thick accretion disk or from the neutron star 
surfac;c, together with hard component which may be 
interpreted as emission from a corona-like structure or 
a boundary layer between the disk and a neutron star 
(White, Stella & Parmar 1988; Christian & Swank 1997; 
Sidoli et al. 2001). We used such two-component models 
to approximate spectra of the three brightest globular clus- 
ter sources showing spectral cut-off (sources ##22, 42, 43). 
and 32) 

For the soft component we used multicolor disk- 
blackbody (XSPEC DISKBB)(Mitsuda et al. 1984) and 
a black body (XSPEC BBODYRAD) models. DISKBB 
model has two parameters, the effective radius, rinVcosi, 
where ri„ is the inner radius of the disk, i is the inclination 
angle of the disk and kTin is the maximum color temper- 
ature of the disk. Although this model does not take into 
account the effects of electron scattering and uses sim- 
plified temperature profile in the accretion disk (Shakura 
& Sunyaev 1973), it has been found to give an adequate 
qualitative description to the spectra of accreting X-ray 
binaries (Mitsuda et al. 1984; White, Stella & Parmar 
1988; Christian & Swank 1997). The BBODYRAD model 
has two parameters, the effective radius, rsB and color 
temperature, kTsB- 

The hard spectral component can be adequately de- 
scribed by various phenomenological and physical models 



involving a break in the slope of the spectrum or quasi- 
exponential spectral cut-off at higher energies. For the 
sake of easier comparison with the results for the Galactic 
LMXB (White, Stella & Parmar 1988; Christian & Swank 
1997; Church & Balucinska-Church 2001), we use a simple 
black body BBODYRAD (in combination with DISKBB 
as a soft component) and a simple power law (in com- 
bination with BBODYRAD as a soft component) to ap- 
proximate hard component in the spectra of brightest GC 
sources in our sample. The results of two-component ap- 
proximation of the spectra of GC X-ray sources are shown 
in Table 7. 

The two-component models give a good fit to the spec- 
tra of the high-luminosity GC X-ray sources in our sam- 
ple. The characteristic temperature of the soft spectral 
component ranges from 0.5 to 1.2 keV, and the average 
contribution of the soft spectral component to the total 
source luminosity is ~ 30%. Interestingly, even in the case 
of brightest M31 GC sources with total isotropic lumi- 
nosities exceeding Eddington limit for the 1.4M0 neutron 
star, the absorption-corrected luminosity of the soft spec- 
tral component never significantly exceeds that limit. In 
conclusion, the two-component approximation of the spec- 
tra of bright M31 GC sources gives spectral parameters 
consistent with high luminosity LMXBs hosting a neutron 
star (White, Stella & Parmar 1988; Christian & Swank 
1997; Church & Balucinska-Church 2001). 

5.1. Spectral distribution of GC X-ray sources 

In order to characterize the overall spectral properties of 
the bright GC X-ray sources, we constructed a distribution 
of their spectral indices in the 0.3 — 7.0 keV energy range 
using the model fits to both XMM-Newton and Chandra 
data with an absorbed simple power law (Fig. 5). For 
the sources with multiple spectral measurements, we used 
average values of the photon index. The main part of the 
distribution including 29 sources out of 31 spans a range 
of photon indices between ~ 1.4 and ^ 2.2, and has a 
clearly defined maximum at a 1.6 (Fig. 5). Excluding 
two extremely hard X-ray sources, gives the corresponding 
weighted mean value of the photon index 1.65±0.01. This 
distribution is reminiscent of the spectral distribution of 
the X-ray sources in the central bulge of M31 (Shirey et al. 
2001). The X-ray population of the central bulge of M31 
is very likely to be dominated by bright LMXBs. Based on 
the similarity in the properties of M31 GC X-ray sources 
and their Galactic counterparts, we expect most of M31 
GC X-ray sources to be LMXBs. Therefore the similarity 
between the spectral distributions of the central bulge of 
M31 and that of its GC X-ray sources is not surprising. 
Moreover, the peak of the M31 GC spectral index dis- 
tribution coincides with the average value of the spectral 
photon index (a ~ 1.6) derived for the LMXB populations 
in the nearby galaxies of different morphological type (Ir- 
win, Athey & Bregman 2003). The spectral photon indices 
of the M31 GC X-ray sources in our sample were found to 
lie in the range, observed for the Galactic GC sources, in 
agreement with previous studies (Di Stefano et al. 2002). 

5.2. Hardness-luminosity diagram of GC X-ray sources 

We studied the relation between the hardness of the 
spectrum and luminosity of GC X-ray sources in our sam- 
ple. In Fig. 6 the hardness of the spectrum of GC X- 
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ray sources expressed in terms of spectral photon index 
is shown as a function of their X-ray luminosity in the 
0.3 — 10 keV energy band. As it is clearly seen from Fig. 
6, there are three groups of sources, occupying distinct re- 
gions in the hardness-luminosity diagram. The first group 
represents majority of the brightest X-ray sources with lu- 
minosities above ~ 4 x 10^^ ergs s~^. The objects in this 
group have relatively narrow range of photon indices with 
average value of 1.55, suggesting common spectral forma- 
tion mechanism at high luminosities. The second group 
comprised of fainter sources, shows greater variety of pho- 
ton indices ranging from ~ 1.4 to ~ 2.3. Based on the wide 
range of X-ray luminosities, this group could represent a 
mix of LMXBs in both low and high luminosity states 
(White, Stella & Parmar 1988). The third group consists 
of two objects with extremely hard spectra (a ~ 0.6 — 0.8). 
Both objects also have a more pronounced spectral cut-off 
than the rest of the GC sources (Table 6). In addition, 
one of the sources, the X-ray source #32 in the globu- 
lar cluster Bo 158, shows regular X-ray dips, implying a 
system observed at high inclination angle (Trudolyubov et 
al. 2002b). The hardness-luminosity diagram of GC X-ray 
sources in our sample appears to be consistent with that 
of the central bulge (Shirey et al. 2001). 

5.3. Low energy absorption 

The unprecedented sensitivity of XMM and Chandra ob- 
servations allows a reliable measurement of low energy ab- 
sorption toward brighter sources in M31. Using the re- 
sults of spectral fitting, we have been able to estimate the 
value of equivalent hydrogen absorbing column, TVh for 31 
source in our sample. To make a homogeneous sample, 
we used the same simple spectral model to estimate A^h 
for all sources (an absorbed power law model). It should 
be noted, that the amount of low-energy absorption de- 
rived from spectral fitting sometimes depends strongly on 
the type of continuum model used to approximate X-ray 
spectrum (Table 4,6). The simple power law approxima- 
tion does not take into account possible presence of an 
additional soft X-ray component leading to the underesti- 
mation of the absorbing column. The choice of a particular 
model for the soft X-ray component (i.e. black body vs. 
multicolor disk black body) or a hard spectral continuum 
(i.e. a Gomptonization model vs. simple power law) in- 
troduces additional uncertainty in the value of N-^. 

Our analysis shows that for the large fraction of GC 
sources the derived value of N^i is either in excess or con- 
sistent with Galactic hydrogen column AT^"' ~ 7 x 10^° 
cm~^ in the direction of M31 (Dickey & Lockman 1990). 
In the first case, it is natural to assume that additional 
absorption could be due to interstellar matter outside our 
Galaxy (i.e. in M31 and intergalactic medium) and/or 
within the X-ray binary itself. We compared the values of 
A'^H, derived from our spectral fits with values predicted 
using the optical observations. If there is no extra ab- 
sorption within the X-ray binary itself, a close agreement 
between the absorbing column values derived from anal- 
ysis of X-ray and optical data is expected (Sidoli et al. 
2001). 

Using the results of optical observations of M31 globular 
clusters (Barmby et al. 2000; Huchra, Brodie & Kent 1991; 
Perrett et al. 2002), we were able to estimate the values of 
optical extinction for 15 GC sources with measured X-ray 



low-energy absorption. The values of optical color excess, 
Eb-v were converted to the equivalent hydrogen column 
densities, N^^ using the relation: Nyi/Eb-v = 5.9 x lO^^ 
cm~^ mag^^ (Cox 2000). In Fig. 7 {left panel) the values 
of the X-ray low-energy absorption, A^^ are plotted against 
optically determined values, A^^''*. The dotted line corre- 
sponds to equal X-ray and optically determined absorp- 
tion. Taking the aforementioned systematic uncertainties 
in the derivation of X-ray absorption into account there 
is a general agreement between two absorption measure- 
ments (Fig. 7, left panel). Based on the current observa- 
tions it is impossible to confirm or rule out a possibility of 
additional absorption in some sources. More sensitive ob- 
servations with wider bandpass (i.e. future Constellation- 
X) and physically justified detailed spectral modeling are 
needed to resolve this interesting presently open issue. 

The distribution of the X-ray absorption could provide 
additional information on the geometry of M31 globular 
cluster system and help to map the structure of ISM inside 
M31. For example, the sources located behind or inside the 
disk of M31 should be in general more absorbed/reddened 
than sources seen in front of the disk. We constructed the 
distribution of absorption columns for 31 M31 GC sources 
(Fig. 7, right panel). The distribution has a prominent 
peak centered at ~ 1.2 x 10^^ cm~^ and a tail structure 
extending up to ^ 6 x 10^^ cm~^. In general, the ob- 
served A^H distribution could be consistent with homoge- 
neous quasi-spherical spatial distribution of GC sources 
with an additional absorption from the disk. 

6. VARIABILITY 

6.1. Long-term variaMlity 

We searched for the long-term flux variability of M31 
GC sources in our sample combining the data of multi- 
ple Chandra and XMM-Newton observations. We found 
that more than 80% of sources with multiple flux mea- 
surements available show significant variability (Table 3). 
In Fig. 8 the long-term flux histories of ten M31 GC 
sources demonstrating the highest levels of variability are 
shown. These sources show both irregular and possible 
quasi-periodic patterns of X-ray variability with corre- 
sponding flux changes by more than a factor of 2 over 
a time scales of several days to years. 

The persistent low luminosity GC X-ray sources with 
average luminosities below few x 10'^'' ergs s"-'^ tend to 
be generally more variable than brighter ones. The cor- 
responding ratio between the highest and the lowest flux 
levels observed with XMM and Chandra is between ~ 3 
and > 15 for low luminosity sources, while it is usually 
lower than ^ 3 for the high luminosity sources. 

The X-ray source associated with globular cluster Bo 
107 shows a possible pattern of recurrent outbursts last- 
ing for ^ 50—100 days spaced by ^ 100 — 200 day intervals. 
The source X-ray flux during the peak of the outburst can 
be at least ~ 10 times higher than during the low lumi- 
nosity periods. This type of variability is reminiscent of 
the Galactic LMXB like Aql X-1 (Priedhorsky & Terrefl 
1984). 

Di Stefano et al. (2002) proposed that the X-ray source 
in Bo 86 may have long-term periodicity on a time scale 
of ~ 200 days (see also Kong et al. 2002). Our data, 
though covering a longer period of time, does not allow us 
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to confirm or completely rule out that kind of variability 
of Bo 86. More frequent observations would be needed to 
establish variability of the source. 

We also searched for the long-term spectral variability 
of all bright M31 GC sources using the results of spec- 
tral analysis of XMM and Chandra data. Several sources 
demonstrate a correlation between the level of X-ray flux 
and hardness of their energy spectrum. The spectral fits 
for most of variable sources indicate that as the source flux 
increases, the spectrum becomes harder (Table 4). Two 
characteristic examples of such behavior are shown in Fig. 
10 (sources Bo 78 and Bo 148). In this figure the hard- 
ness of X-ray spectrum expressed in terms of best-fit pho- 
ton index is shown as a function of the source luminosity. 
The luminosity of the spectrally variable sources ranges 
from ~ 10^^ to ~ few x 10"^® ergs s~^, falling into typical 
range of luminosities for both atoll and Z-source classes 
of Galactic LMXB (Hasinger & van der Klis 1989). For 
brighter sources with luminosities around 10^^ ergs s~^ the 
observed spectral changes could be consistent with that of 
the Z-sources on the normal branch (Hasinger & van der 
Klis 1989). 

6.2. RecAirrent transient activity of Bo 163. 

One of the sources, coincident with a globular cluster Bo 
163 shows a recurrent transient-like variability. Our anal- 
ysis of the archival Chandra data revealed the presence 
of a bright X-ray source CXO J004317.8+412745 with a 
position consistent with Bo 163 and previously detected 
ROSAT source RX J0043.2-I-4127 (Supper et 'al. 2001) 
during 2001 February 18 ACIS-S observation (Obs. ID 
#1582) (Fig. 9, upper panel). The estimated X-ray lumi- 
nosity of the source was ^ 10"^® ergs s^^ in the 0.3 — 10 
keV energy band (see source spectral parameters in Table 
5). Other XMM-Newton and Chandra observations cov- 
ering the position of Bo 163 did not yield positive source 
detection with 2a upper limits on its luminosity ranging 
from ~ 10"^^ to ~ 10^^ ergs s~^ (depending on the type of 
instrument and exposure time). 

We analyzed the data of the earlier observations of Bo 
163/RX J0043. 2+4127 with ROSAT/PSPG and HRI de- 
tectors. The observed ROSAT source count rates were 
converted into the X-ray fluxes in the 0.3—10 keV energy 
band using WebPIMMS assuming an absorbed power law 
spectrum with iVn = 7 x 10^" cm^^ and a = 1.7. The re- 
sulting ROSAT light curve of RX J0043. 2+4127 is shown 
in lower panel of Fig. 9. The source clearly demonstrates 
a pattern of recurrent transient-like behavior with bright 
outbursts with Lx > 10^^ ergs and an extended qui- 
escent periods with Lx < 10'^® ergs s~^. The correspond- 
ing maximum ratio of the the highest outburst luminosity 
(Lx ~ 2 X 10"^^ ergs s~^) to the lowest measured quiescent 
luminosity {Lx < 10^^ ergs s"^) is > 2000. 

The overall X-ray properties of RX J0043.2+4127 are 
reminiscent of both Galactic neutron star and black hole 
recurrent transients like Cen X-4 (Kaluzienski, Holt & 
Swank 1980, van Paradijs et al. 1987) and 4U1630-40 
(Priedhorsky 1986). Unfortunately, current spectral and 
timing data does not allow us to determine the nature 
of the compact object in RX J0043. 2+4127. Observations 
with higher sensitivity and timing resolution are needed to 
choose between neutron star and black hole interpretation. 



6.3. Short-term variability. 

The comparisons based on the broad-band spectral 
properties are not sufficient to establish a neutron star 
nature or rule out a black hole nature for the M31 globu- 
lar cluster sources. On the other hand, the study of their 
short-term variability can provide a definitive answer, if 
Type I X-ray bursts or X-ray pulsations are observed. We 
searched for both types of variability in the M31 GC data. 
Unfortunately, the observed source count rates for most 
sources are too low to allow us to probe timescales shorter 
than ~ 20 — 30 seconds. We did not find evidence of X- 
ray pulsations and short X-ray bursts in the XMM-Newton 
and Chandra data. 

Several GC sources demonstrate significant aperiodic 
variability on a time scale of individual observation. As 
an example of such behavior, the lightcurves of three GC 
sources arc shown in Fig. 11. All three sources (Bo 86, Bo 
107 and Bo 135) show irregular changes of X-ray fiux of 
up to ~ 100% over a time scales of 1000 - 10000 s. These 
sources were also found to be variable on a time scales of 
months to years (Fig. 8). The combination of such short 
and long term variabilities was long recognized as one of 
the characteristic features of X-ray binaries (Hasinger & 
van der Klis 1989). In addition, the X-ray source coin- 
cident with M31 GC candidate Bo 158 shows a remark- 
able dip-like modulation with a a period of 10017 s (Tru- 
dolyubov et al. 2002b). 

The unprecedented sensitivity of XMM-Newton provides 
a unique opportunity to study spectral evolution of the 
brightest LMXB sources in M31 on a time scale of hun- 
dreds and thousands of seconds. To study short-term spec- 
tral variability of the bright GC sources, we constructed 
their X-ray hardness-intensity diagrams using the data of 
XMM-Newton observations. The X-ray hardness was de- 
fined as the ratio of the source intensities in the 2.0 — 7.0 
keV and 0.3 — 2 keV energy bands with data integration 
times of 1000 — 3000 s depending on the source intensity. 
Several sources show a complex patterns of evolution on 
the hardness-intensity diagram somewhat reminiscent of 
the Galactic Z and atoll sources (Hasinger & van der Klis 
1989). In Fig. 12, we show the hardness-intensity dia- 
gram of one of these sources, the X-ray source in globular 
cluster Bo 148. Although the data strongly suggest some 
pattern of spectral evolution, it is still insufficient to clas- 
sify source behavior in conventional framework of Galactic 
LMXB. More high-sensitivity observations covering wider 
range of source luminosities are needed to address this is- 
sue. 

6.4. "Dipping" X-ray source in Bo 158 

The X-ray source J004314.1+410724 was discovered in 
M31 by the Einstein observatory (source #85 in Trinchieri 
& Fabbiano 1991) and was detected in subsequent ob- 
servations with ROSAT (Primini et al. 1993, Supper et 
al. 2001), XMM-Newton (Trudolyubov et al. 2002b) and 
Chandra (Di Stefano et al. 2002). The position of the 
source is consistent with the optically identified globu- 
lar cluster candidate Bo 158 (source #158 in Table IV 
of Battistini et al. 1987). The observation of the cen- 
tral bulge of M31 with XMM-Newton taken 2002 Jan 6 
revealed strong periodic dip-like modulation of the source 
flux (Trudolyubov et al. 2002b; Fig. 13, upper panel). The 
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X-ray flux was modulated by ^ 83% at a period of 2.78 
hr (10017 s). The X-ray intensity dips show no energy 
dependence. The analysis of earlier archival XMM and 
i? 05^ T observations taken in the year 2000 and 1991 re- 
vealed weaker dips with amplitudes of ^ 30% and ^ 50%. 
The amplitude of the modulation has been found to be 
anticorrelated with source X-ray flux. A detailed descrip- 
tion of the timing and spectral properties of Bo 158 are 
presented in Trudolyubov et al. 20026. 

We extended our original study of the X-ray modula- 
tion in Bo 158 combining the data of archival XMM and 
Chandra observations. Because of a relatively long period 
of modulation (~ 10000 s), only one Chandra observation 
was selected for our analysis (Obs. #2017 with ACTS). 
For one XMM observation (Obs. #1), showing signiflcant 
change of the source flux on a time scale of ~ 5 hours, 
we divided the data into two segments corresponding to 
different luminosity levels. A sample light curves of Bo 
158 obtained with XMM and Chandra are shown in the 
upper and middle panels of Fig. 13. The modulation frac- 
tion (i.e. the ratio of the average modulation amplitude to 
the average source flux) was determined for each observa- 
tion. The resulting dependence of the modulation fraction 
on the source flux is shown in the bottom panel of Fig. 

13. The relative amplitude of X-ray modulation declines 
dramatically (at least by factor of ~ 10) as the source lu- 
minosity rises from ~ 4.5 x 10^^ to ~ 1.4 x 10^* ergs s~^. 
The observed correlation between the modulation fraction 
and X-ray flux of Bo 158 allows us to put constraint on 
possible number of bright X-ray sources presently active 
in this cluster. We conclude that one X-ray source must 
contribute up to 99% of the detected flux. 

In addition to the short-term periodic variability, there 
is clear evidence for signiflcant long-term variability of Bo 
158. The long-term monitoring observations with Chandra 
and XMM-Newton (Fig. 8) show that the source is highly 
variable on a time scales of months to years. The X-ray 
luminosity Bo 158 changes at least by factor of ~ 4, e.g. 
between ~ 5 x 10^^ and ~ 2 x 10^* ergs s~^. 

7. LUMINOSITY DISTRIBUTION OF M31 GO X-RAY 
SOURCES. 

We have built differential and cumulative luminosity dis- 
tributions of M31 GC X-ray sources assuming a distance 
of 760 kpc (Fig. 14). The luminosity distribution for the 
central bulge based on the results of XMM-Newton obser- 
vations (Trudolyubov et al. 2002 a) is also shown in Fig. 

14. Luminosities of the brighter sources (Lx > 5 x 10^^ 
ergs s~^) were derived using analytical flts to their spectra 
in the 0.3 — 10 keV energy band. For faint sources, lumi- 
nosities were derived using WebPIMMS, assuming an ab- 
sorbed power law model with Nu = 7 x 10^° cm~^ (Dickey 
& Lockman 1990) and a = 1.7 (Shirey et al. 2001) in the 
same energy band. For sources with multiple flux mea- 
surements available, we used average flux to calculate the 
luminosity. Taking into account the sensitivity as a func- 
tion of off-axis distance, different exposure times and the 
effect of diffuse X-ray emission near the center of M31, we 
estimate flux completeness limit of our samples as ^ 10^^ 
ergs s~^ (indicated by dotted line in Fig. 14). 

The shapes of both cumulative and differential luminos- 
ity distributions for M31 GC sources indicate the presence 
of a cut-off and can be described by a broken power law or 



a c\it-off power law models. We fitted the cumiilative lu- 
minosity distribution with two power laws at hmiinositic;s 
below and above 5 x 10'^'' ergs s~^. For source luminosities 
between 10^^ and 5 x 10^^ ergs s~^, we obtain the slope of 
—0.3 ± 0.1 with normalization of 38 sources at 10^^ ergs 
s^^, for higher luminosities, the integral slope is — 1.2^q7 
(shown as dotted lines in Fig. 14). We also fitted the un- 
binned differential luminosity with a power law with expo- 
nential cut-off: i~"exp(— L/Lcut)- The resulting values of 
a and the cut-off luminosity are 1.0 and 2.7l:tl x 10^^ 

ergs 

s~^ respectively. There is some flattening of the cumula- 
tive distribution toward the faint luminosity end, probably 

caused by incompleteness of our sample. 

As it is clearly seen from Fig. 14, the integrated lu- 
minosity function of GC sources differs significantly from 
that of the other point sources in the central bulge of M31. 
The low luminosity part of the GC luminosity distribution 
is somewhat flatter than that of the central bulge of M31. 
The luminosity of the break in the slope of the luminos- 
ity function appears to be ~ 2 — 3 times higher for GC 
sources {L^^ ^ 5 x lO''^ ergs s~^) than for the sources in 

the central bulge with L^^'^^" - (1.5 - 2.0) x lO^^ ergs s'^ 
(Shirey et al. 2001, Trudolyubov et al. 2002a, Kong et al. 
2002 a). Qualitatively, the fraction of the bright sources 
{Lx > few X 10^'' ergs s~^) is higher for the GC popula- 
tion than for the central bulge population. It is interesting 
to note that in spite of the aforementioned differences, the 
slopes of the high-luminosity parts of the luminosity dis- 
tributions for the M31 GC and central bulge populations 
are very close to each other. 

In addition, there is an obvious difference between the 
luminosity functions of GCs and X-ray sources in the 
disk of M31 (Trudolyubov et al. 2002a). The luminos- 
ity distribution of the disk X-ray source population ter- 
minates at much lower luminosities {few xlO^^ ergs s~^) 
(Trudolyubov et al. 2002a) than that of the M31 GCs 
{several x 10'^* ergs s~^). 

It is also interesting to compare the luminosity functions 
of the Galactic GC and M31 GC X-ray sources. These two 
functions have two main differences: i) the peak X-ray lu- 
minosity of the M31 GC X-ray sources is ~ 10 times higher 
than that of the brightest Galactic GC X-ray source (Di 
Stefano et al. 2002); ii) the fraction of the bright sources 
is much higher for the M31 GC population than for the 
Galactic GC population; 13 sources or ~ 30% of the total 
number of GC sources detected in our survey reached lu- 
minosity of 5 x 10''^ ergs s^^, while only one Galactic GC 
source, X1820-303 (NGC 6624) has luminosity occasion- 
ally exceeding this level (Sidoli et al. 2001). 

8. ON THE ROLE OF THE BRIGHT GC SOURCES IN THE 
X-RAY SOURCE POPULATION OF M31. 

We found that ^ 10% of the optically identifled GC can- 
didates in the regions of M31, covered by XMM-Newton 
and Chandra observations, harbor bright X-ray sources 
with luminosities above 10^^ ergs s~^ (43 X-ray sources 
were found among ~ 400 GC candidates). This fraction 
is comparable to the observed for the Milky Way glob- 
ular cluster system with 14 bright X-ray sources among 
known 150 globular clusters (Verbunt et al. 1995, Sidoli 
et al. 2001, Harris 1996). In addition, the fraction of M31 
GCs that contain X-ray sources with luminosities above 
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~ 10^^ ergs 5%) appears to be consistent with that 

derived for a large sample of galaxies of various morpho- 
logical types (Sarazin et al. 2003). 

Our observations clearly indicate that GC X-ray sources 
make dominant contribution to the bright {Lx > few x 
2^q37 (,j.gg g-i^ source counts in the surveyed areas of 
M31. Six persistent globular cluster X-ray sources (sources 
##2,9,22,41,42,43) have absorbed luminosities in the 

0. 3 — 10 keV energy band exceeding 10^^ ergs s~^; 
three other persistent GC sources in our sample (sources 
##31,32,37) occasionally exceed that luminosity level, 

1. e. at least six GC X-ray sources are brighter than 10^* 
ergs s~^ at any given time. Apart from transient X-ray 
sources, we detected only one persistent non-GC source 
belonging to M31; having average luminosity in excess of 
10^8 ergs s-i, RX J0042.6+4115 (Barnard^ Kolb & Os- 
borne 2003), and two other non-GC sources occasionally 
exceeding that luminosity level in the same regions. There- 
fore, from ~ 67% to ~ 90% of the sources brighter than 
2()-5s ^jj,gg g-i observed in the surveyed area at any given 
time are associated with globular clusters. According to 
the results of our survey, about 20 sources or ~ 40% of the 
total number of sources with luminosities above lO''^ ergs 
s~^ detected in XMM and Chandra observations of M31, 
reside in globular clusters. 

We found that fraction of bright X-ray sources {Lx > 
10^'' ergs s~^), residing in globular clusters is signifi- 
cantly higher in M31 (~ 40%) than in Milky Way Galaxy 
(~ 5 - 10%) (Liu et al. 2001). It should be also noted, 
that no persistent GC sources with luminosities in excess 
of 10^* ergs s~^ have been detected in our Galaxy so far, 
compared to six sources found in M31. On the other hand, 
the fraction of bright GC X-ray sources in M31 is surpris- 
ingly similar to the observed in early-type galaxies: i.e. 
NGC 1399 (Angelini, Loewenstein & Mushotzky 2001) and 
NGC 4472 (Kundu, Maccarone, & Zepf 2002). 

9. SPATIAL DISTRIBUTION AND LUMINOSITY OF M31 GC 
X-RAY SOURCES. 

The spatial distribution of GC X-ray sources detected 

with XMM-Netwon and Chandra is shown in Fig. 15. We 
have noticed that bright X-ray sources tend to reside at sig- 
nificantly larger galactocentric distances than the fainter 
sources. In fact, all GC X-ray sources with luminosities 
steadily exceeding 10^^ ergs s~^ (shown in Fig. 15 with 
large white circles) happen to have angular distances more 
than 15' from the center of the galaxy, or well outside the 
central bulge of M31, where the most of fainter GC X-ray 
sources are residing. 

10. ON THE CORRELATION BETWEEN THE X-RAY AND 
OPTICAL PROPERTIES. 

We used the results of optical observations of M31 GC 

candidates (Huchra, Brodie & Kent 1991; Barmby et al. 
2000; Perrett et al. 2002) to study how the properties of 
the X-ray sources depend on the optical properties of the 
globular clusters hosting them. 

We found that globular clusters hosting bright X-ray 
sources tend to be optically brighter than the rest of M31 
globular clusters, in agreement with previous studies (Di 
Stefano et al. 2002). To illustrate this effect, V-magnitude 
distribution of X-ray bright globular clusters from our 
sample is shown in Fig. 16{left panel) along with the 



distribution for the whole sample of GC candidates from 
Battistini et al. (1987). Similar tendency for bright X- 
ray sources to be associated with more optically luminous 
GCs has been also observed in several early-type galaxies 
(Sarazin et al. 2003). Sarazin et al. (2003) shows that 
this tendency could be consistent with a constant proba- 
bility of finding an LMXB per unit optical luminosity of 
the cluster. 

Both Galactic and M31 globular clusters hosting bright 
LMXBs were found to be both denser and more metal-rich 
(Bellazzini et al. 1995). We studied the effect of globular 
cluster metallicity, [Fe/H] on the properties of the M31 
GC X-ray sources from our sample. For 31 GC X-ray 
sources the results of metallicity measurements are avail- 
able (Huchra, Brodie & Kent 1991; Barmby et al. 2000; 
Perrett et al. 2002). The metallicity distribution of the 
X-ray bright GC in our sample is shown in Fig. 16{right 
panel). The population of X-ray emitting GC yields a 
weighted mean < [Fe/H] >xr= — 1.06±0.02, significantly 
higher than < [Fe/H] >nXR= —1.37 ±0.01 derived for the 
non- X-ray GC subsample from Perrett et al. (2002), in 
agreement with results of previous studies (Bellazzini et 
al. 1995). We used a Kolmogorov-Smirnov (KS) test to 
determine whether X-ray emitting GC and non- X-ray GC 
populations were drawn from the same distribution. The 
hypothesis that the two distributions are extracted from 
the same parent population can be rejected at the confi- 
dence level of ^ 93%. 

The dependence of the X-ray Imninosity of M31 GC 
sources on the metallicity of their host globular clusters 
is shown in the left panel of Fig. 17. The data points 
can be separated into two main groups corresponding to 
the source X-ray luminosities below and above ~ 10'^* 
ergs s~^. The first group (sources with X-ray luminosi- 
ties below 10^^ ergs s~^) has weighted mean metallicity of 
—0.79 ± 0.04 and possibly demonstrates general tendency 
of the brighter sources to reside in the more metal rich 
clusters. This behavior could be consistent with possible 
correlation between X-ray luminosity and cluster metallic- 
ity found by Sidoli et al. (2001) for the sample of bright 
Galactic GC sources. It is also interesting to note, that 
all known Galactic GC X-ray sources have luminosities in 
the range occupied by the first group (Sidoli et al. 2001). 
The whole second group (sources with X-ray luminosities 
above 10^^ ergs s~^) occupies low metallicity region with 
weighted mean metallicity of —1.12 ± 0.02. 

Recent X-ray observations and theoretical work suggest 
that spectral hardness of the bright LMXB sources in the 
soft X-ray band (0.1 — 2.4 keV) can depend on the average 
metallicity of the globular clusters harboring them. The 
soft X-ray spectra of the sources located in the metal-poor 
GCs were found to be harder than those of the metal-rich 
clusters (Irwin & Bregman 1999; Maccarone, Kundu & 
Zepf 2004). To demonstrate the effect of cluster metal- 
licity on the intrinsic spectral hardness of the M31 GC 
X-ray sources, in the right panel of Fig. 17 the spectral 
power law photon index in the 0.3 — 7 keV energy band for 
the GC sources from our sample is shown as a function of 
metallicity of the globular clusters hosting them. As it is 
seen from Fig. 17 [right panel), there is a possible corre- 
lation of the X-ray spectral hardness and the metallicity. 
We note however, that it is premature to draw any con- 
clusions on the relationship between the intrinsic spectral 
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parameters of X-ray sources and metallicity of their host 
clusters, based on a model-dependent presentation of the 
spectrum measured in the narrow energy band. For ex- 
ample, using a multi-component spectral modeling of the 
BeppoSAX broad-band data for the bright Galactic GC 
sources, Sidoli et al. (2001) foimd a correlation between 
the tcmperatiirc of the soft spectral component and a clus- 
ter metallicity, opposite to the results obtained in the soft 
X-ray band with one-component spectral models (Irwin & 
Bregman 1999). 

11. ON THE NATURE OF M31 GC X-RAY SOURCES 

Several GC X-ray sources observed in M31 have lumi- 
nosities exceeding 10"^* ergs s~^, being much brighter than 
any of the Galactic GC X-ray sources (Table 3). Four 
X-ray sources in our sample (##2,10,22,42) have lumi- 
nosities close to or exceeding a formal Eddington limit for 
the 1.4Mq neutron star. Nevertheless, such a high lu- 
minosities are not uncommon for bright low-mass X-ray 
binaries containing neutron star primaries. The Z-source 
subclass of Galactic LMXB is characterized by high X- 
ray luminosities (> 10^^ ergs s~^) (Hasinger fc van der 
Klis 1989) with two sources, Sco X-1 and GX 5-1, regu- 
larly exceeding Eddington limit for a neutron star (Liu et 
al. 2001). Recently, Barnard, Kolb and Osborne (2003) 
found that bright M31 X-ray source, RX J0042.6-h4115 
with extremely high X-ray luminosity of ~ 5 x 10^* ergs 
could be also classified as a neutron star Z-source. On the 
other hand, among stellar mass black hole X-ray sources, 
only transients in outburst (Chen, Shrader & Livio 1997) 
and highly variable systems like GRS 1915-1-105 (Morgan, 
Remillard & Greiner 1997) reach luminosities observed for 
bright M31 GC sources. 

All bright X-ray sources found in the Milky Way glob- 
ular clusters show Type I X-ray bursts and probably con- 
tain neutron star primaries. Moreover, theoretical cal- 
culations predict the probability of finding a black hole 
X-ray binary in the globular cluster to be much smaller 
than the neutron star binary due to dynamical effects and 
extremely low duty cycles of black hole binaries (Kulka- 
rni et al. 1993, Sigurdsson & Hernquist 1993, Kalogcra, 
King & Rasio 2004) . The considerations mentioned above, 
and a general similarity between the M31 and Milky Way 
globular cluster systems allow us to expect most of the 
persistent M31 GC X-ray sources to be LMXB systems 
containing a neutron star primaries. However, the possi- 
bility of finding an active accreting black holes in our GC 
source sample can not be ruled out, based on the current 
spectral and timing data. For example, the X-ray source 
RX J0043.2-h4127 associated with globular cluster Bo 163 
shows a pattern of recurrent transient outbursts and qui- 
escent periods (Fig. 9) typical for both neutron star and 
black hole transient systems. The study of short-term vari- 
ability of GC X-ray sources in M31 with future missions 
can provide a definitive answer, if Type I X-ray bursts 
or X-ray pulsations are observed from these objects. The 
broad-band long-term monitoring of the transient /highly 
variable GC sources could also provide valuable informa- 
tion on their spectral evolution helping to establish the 
nature of these systems. 

Some of the M31 GC X-ray sources in our sample might 
be unresolved composites of separate sources (Angelini, 
Loewenstein & Mushotzky 2001; Di Stefano et al. 2002). 



Recent high-resolution Chandra observations of the Galac- 
tic globular clusters revealed that one of the bright Galac- 
tic GC sources is a composite of two X-ray sources of com- 
parable brightness (White & Angelini 2001). Depending 
on the probability of finding a bright X-ray source in a GC, 
we can expect to find 2 — 4 GC X-ray sources in our sample 
(or ~ 5 — 10% of the total number) to be composites of 
two independent sources (Di Stefano et al. 2002). 

12. SUMMARY 

Utilizing the data of XMM-Newton and Chandra ob- 
servations, we performed a study of individual and group 
properties of 43 X-ray sources identified with globular clus- 
ters in M31. The correlation between the optical param- 
eters of the host globular clusters and the properties of 
X-ray sources associated with them was also studied. We 
compare the individual and group properties of M31 GC 
X-ray sources with the properties of their Galactic coun- 
terparts. 

The spectral properties of the bright GC X-ray sources 
in our sample were found to be similar to that of the 
LMXB located in the bulge and globular clusters of the 
Milky Way Galaxy at the same luminosity levels (White, 
Stella & Parmar 1988; Christian & Swank 1997; Church 
& Balucinska-Church 2001). We note, that some GC 
sources could be unresolved composites of two or even 
more sources of comparable brightness complicating direct 
comparison with their Galactic counterparts. For the ma- 
jority of sources, we obtained acceptable fits using a power 
law spectral model approximation. However, for several 
sources with high luminosities, the models with quasi- 
exponential cut-off at ^ 2 — 8 keV or two-component mod- 
els describe the energy spectra significantly better than a 
simple power law. Several sources demonstrate a correla- 
tion between the level of X-ray flux and hardness of their 
energy spectrum similar to the Galactic Z and atoll-sources 
(Hasinger & van der Klis 1989). The spectral distribution 
of M31 GC X-ray sources resembles the corresponding dis- 
tributions derived for the the central bulge of M31 (Shirey 
et al. 2001) and other nearby galaxies of different mor- 
phological type (Irwin, Athey & Bregman 2003). 

There is no evidence that X-ray variability of M31 GC 
sources and bright Galactic LMXB (both GC and non-GC) 
differ. We found that ~ 80% of the sources in our sample 
with multiple flux measurements available show significant 
variability on a time scales from days to years. We note, 
that in some GC X-ray sources the observed variability 
could be result of superposition of flux variations of several 
unresolved sources. The X-ray source RX J0043. 2-1-4127 
in Bo 163 has been found to show recurrent transient out- 
bursts with peak luminosities above lO'"*^ ergs s^"'^, typical 
for Galactic recurrent transients that harbor both neutron 
stars and black holes. Several sources show significant vari- 
ability on a time scale of individual observations, ranging 
from aperiodic fluctuations to regular dipping (Bo 158). 

The X-ray luminosity function of GC sources is found 
to be significantly different from that of the other point 
sources in the bulge and disk of M31 and Galactic GC 
X-ray sources. The luminosity distributions of the GC X- 
ray sources in M31 and that of the Milky Way have two 
main differences: i) the peak luminosity of the M31 GC 
X-ray sources is 10 times higher than that of the bright- 
est Galactic GC X-ray source; ii) the fraction of the bright 
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sources is much higher for the M31 GC population; ^ 30% 
of the GC sources detected in our survey reached the hi- 
minosity of 5 x 10^^ ergs s~^, while only one of 14 Galactic 
GC sources, X1820-303 (NGC 6624) has luminosity occa- 
sionally exceeding this level (Sidoli et al. 2001). 

We found that ~ 10% of the optically identified GC can- 
didates in the regions of M31, covered by XMM-Newton 
and Chandra observations, harbor bright X-ray sources 
with luminosities above 10^^ ergs s~^. This fraction is 
comparable to that observed for the Milky Way globular 
cluster system with 14 bright X-ray sources among known 
~ 150 globular clusters (Verbunt et al. 1995, Sidoli et 
al. 2001, Harris 1996). In addition, the fraction of M31 
GCs that contain X-ray sources with luminosities above 
^ 10'^''' ergs s^^ 5%) appears to be consistent with that 
derived for a large sample of galaxies of different morpho- 
logical type (Sarazin et al. 2003). 

The inferred isotropic X-ray luminosities of the GC 
sources lie between ~ 10^^ and ~ 10^^ erg s~^ in the 
0.3 — 10 keV energy band. Six persistent sources (or ~ 14% 
of the total number) have luminosities steadily exceeding 
10^^ erg s~^, and three other sources occasionally exceed 
that luminosity level. Four X-ray sources in our sample 
(##2,10,22,42) have luminosities close to or exceeding 
the Eddington limit for the 1.4M0 neutron star. Our ob- 
servations indicate that GC sources make dominant contri- 
bution to the bright source counts in the areas of M31 cov- 
ered by the survey: ^ 40% of the total number of sources 
with luminosities above 10"^^ ergs s^^ reside in GCs. The 
fraction of the GC X-ray sources rises to 67 — 90% if one 
considers source luminosities above 10^^ erg s~^. The con- 
tribution of the GC sources to the total number of bright 
sources found in M31 is much higher (by factor of > 4) 
than in the Milky Way galaxy. On the other hand, there 
is a surprising similarity between the fractions of bright 
GC X-ray sources in M31 and in the early-type galaxies 
(Angelini, Loewenstein & Mushotzky 2001, Kundu, Mac- 
carone, & Zepf 2002). 

The distribution of radial distances of GC X-ray sources 
from the center of M31 shows that brightest systems 
with luminosities steadily exceeding 10^^ erg s"""^ reside 
at greater galactocentric distances than the rest of the GC 
sources. 

We used the results of optical observations of M31 
GC candidates to study how the properties of the X-ray 
sources depend on the optical properties of the globular 



clusters hosting them. We found that globular clusters 
hosting bright X-ray sources tend to be optically brighter 
and more metal rich than the rest of M31 globular clus- 
ters, in agreement with previous studies (Bellazzini et al. 
1995, Di Stefano et al. 2002). The comparison of the X- 
ray properties of GC sources with metallicity of their host 
globular clusters has shown that the brightest sources with 
luminosities above ~ 10^^ ergs s~^ tend to reside in more 
metal poor clusters than the fainter sources. 

The remarkable similarities between the properties of 
the M31 GC X-ray sources and that of the Galactic neu- 
tron star LMXBs allow us expect most of the persistent 
M31 GC X-ray sources to be LMXB systems with neu- 
tron star primaries. We note, however, that current X-ray 
spectral and timing data can not rule out the possibility 
of finding an active accreting black holes in our GC source 
sample. The study of short-term variability and broad- 
band spectral evolution of GC X-ray sources in M31 with 
future missions can provide valuable information allowing 
to establish the nature of these systems. 

In conclusion, we note that in spite of the similarities 
of the spectral shape and variability patterns of individ- 
ual M31 GC X-ray sources and Galactic LMXBs (both 
GC and non-GC), the group properties of M31 and Milky 
Way globular cluster X-ray sources differ significantly. The 
M31 GC source population has a much higher fraction of 
bright X-ray sources and higher peak luminosity than the 
Milky Way GC population, probably suggesting different 
formation histories and evolution. The large fraction of 
GC systems among the bright X-ray sources in M31 makes 
the M31 population of GC X-ray sources more like giant 
early-type galaxies than the Milky Way. 
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Table 1 

XMM-Newton observations of M31 fields used in our analysis. 



Obs. # 


Date, UT 


Field 


Obs. ID 


RA (J2000)" 


Dec (J2000)" 


Exp.(MOS)'' 


Exp.Cpn)" 










(h:m:s) 


(d:m:s) 


(ks) 


(ks) 


#1 


2000 Jun 25 


M31 Core 


0112570401 


00:42:43.0 


41:15:46.1 


28.9 


24.9 


#2 


2000 Dec 28 


M31 Core 


0112570601 


00:42:43.0 


41:15:46.1 


12.1 


9.4 


#3 


2001 Jan 11 


M31 CI 


0065770101 


00:32:46.9 


39:34:41.7 


7.3 


4.9 


#4 


2001 Jun 29 


M31 Core 


0109270101 


00:42:43.0 


41:15:46.1 


29.0 


24.9 


#5 


2002 Jan 05 


M31 Northl 


0109270701 


00:44:01.0 


41:35:57.0 


57.3 


54.7 


#6 


2002 Jan 06 


M31 Core 


0112570101 


00:42:43.0 


41:15:46.1 


63.0 


49.9 


#7 


2002 Jan 13 


M31 Southl 


0112570201 


00:41:25.0 


40:55:35.0 


53.9 


53.3 


#8 


2002 Jan 24 


M31 South2 


0112570301 


00:40:06.0 


40:35:24.0 


30.0 


29.0 


#9 


2002 Jan 26 


M31 North2 


0109270301 


00:45:20.0 


41:56:09.0 


29.1 


25.3 


#10 


2002 Jun 29 


M31 North3 


0109270401 


00:46:38.0 


42:16:20.0 


51.5 


36.5 



" - coordinates of the center of the field of view 



'' - instrument exposure used in the analysis 
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Table 2 

Chandra observations of M31 fields used in our analysis. 
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1584/ACIS-I 


00:42:37.50 
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4.9 


2001-07-03 


M31-3 


2054/ACIS-S 
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41:40:55.0 


14.7 


2001-07-03 


M31-2 
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00:41:49.90 


40:59:20.0 


13.8 
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00:42:42.30 
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46.9 


2001-11-19 


M31 NUCLEUS 
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00:43:05.55 


41:17:03.3 
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00:43:09.80 


41:19:00.7 


4.9 


2002-02-06 


M31 TRANSIENT 


2896/ACIS-I 


00:43:05.50 


41:17:03.3 


4.9 


2002-06-02 


M31 TRANSIENT 


2898/ACIS-I 


00:43:09.80 


41:19:00.7 


4.9 


2002-07-08 


M31 TRANSIENT 


2901/ACIS-I 


00:41:54.67 


40:56:47.5 


4.6 


2002-08-11 


M31 


4360/ACIS-I 


00:42:44.40 


41:16:08.9 


4.9 


2002-08-23 


M31 TRANSIENT 


2899/ACIS-I 


00:41:54.69 


40:56:47.8 


4.9 


2002-10-14 


M31 TRANSIENT 


2894/ACIS-I 


00:42:34.90 


40:57:21.0 


4.7 



" - coordinates of the center of the field of view 



* - total exposure 
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Table 3 

List of GC X-ray sources detected in XMM-Newton and Chandra observations of M31. 



Source 
ID 




R.A. 


Dccl. 


Lx(0.3-10 keV) 
(xlO^'^ ergs s"^) 


Optical ID" 


X-ray ID" 


1 


00 






oy 


^4 
o^ 


40 


6 


Gl 




2 


00 








^o 


O I 


1990 


Bo5 


SHP73 


3 


00 


41 


41.0 


41 


04 


01 

U J. 


3 


MIT87 


D27 


4 


00 


41 




41 


12 


12 


2.4-2.9 


Bo55 




5 


00 


41 


^19 Q 


An 


47 


1 n 


1-24 


Bo58,MIT106 


D22 


6 


00 


42 




41 


09 
uz 


48 


31-73 


Bo D42,MIT130 


SHP138,D13 
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00 




07 1 


41 


00 


1 fi 
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00 
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1 7 
J. I 


4'i 


40-114 


MIT 140 


DIO 


9 


00 


49 


19 1 
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1 7 


oo 
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01 
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00 
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41 


14 


01 
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12 
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42 


1 Q (S 


41 


21 


oo 
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00 


42 
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"17 


1 Q 


19 
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14 


00 


49 
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41 


1 Q 
±y 


1 ^ 

lO 


28-174 


Bo96,MIT174 


D14 


15 


00 


49 


97 4 


An 


oy 


ou 


1-12 


Bo98 


D25 


16 


00 


42 


31.2 


41 


1 Q 

J- y 


■^8 

oo 


56-290 


Bol07,MIT192 


SHP175,D16 


17 


00 


42 


oo. J. 


41 


0^ 
uo 


98 
zo 


41-52 


BollO 


SHP178 


18 


00 


42 


34.4 




^7 

O I 


OQ 
uy 


14 


Boll7 




19 


00 


42 




41 


10 


^9 
oz 


16-27 


Bol23,MIT212 


D18 


20 


00 


42 


41.4 


41 


1 ^ 


9^ 
zo 


34-137 


MIT213 


D23 


21 


00 


49 


"ifl 7 


41 


1 


oo 


3 


MIT222 


D28 


22 


00 


49 


'il Q 
o ± . y 


41 


■^1 


07 


3093-4009 


Bol35 


SHP205 


23 


00 


49 


oo.o 


41 


1 8 


oo 


8-84 


Bol38 




24 


00 


49 


oy.o 


A 1 


ly 


1 Q 

ly 


152-555 


Bol43 


SHP217,D5 


25 


00 


42 


59.8 


41 


16 


05 


216-512 


Bol44 


D6 


26 


00 


43 


01.3 


41 


30 


17 


553-692 


Bo91,MIT236 


SHP218 


27 


00 


43 


02.9 


41 


15 


22 


74-414 


Bol46 


SHP220,D7 


28 


00 


43 


03.3 


41 


21 


22 


52-194 


Bol47,MIT240 


SHP222,D11 


29 


00 


43 


03.8 


41 


18 


05 


105-418 


Bol48 


SHP223,D8 


30 


00 


43 


07.5 


41 


20 


20 


16-72 


Bol50,MIT246 


D19 


31 


00 


43 


10.6 


41 


14 


50 


373-1248 


Bol53,MIT251 


SHP228,D3 


32 


00 


43 


14.2 


41 


07 


24 


600-1880 


Bol58 


SHP229,D12 


33 


00 


43 


14.7 


41 


25 


13 


2 


Bol59 




34 


00 


43 


15.4 


41 


11 


24 


16-22 


Bol61,MIT260 


D26 


35 


00 


43 


17.8 


41 


27 


45 


< 1 - 1010 


Bol63 


RX J0043.2-F4127 


36 


00 


43 


36.7 


41 


08 


10 


46 


Bol82 




37 


00 


43 


37.3 


41 


14 


43 


454-1981 


Bol85,MIT299 


SHP247,D9 


38 


00 


43 


42.9 


41 


28 


50 


52 


MIT311 


SHP250 


39 


00 


43 


45.5 


41 


36 


57 


44 


Bol93 


SHP253 


40 


00 


43 


56.5 


41 


22 


02 


37 


Bo204 


SHP261 


41 


00 


44 


29.6 


41 


21 


36 


1130 


Bo225 


SHP282 


42 


00 


45 


15.4 


11 


39 


42 


5148-10372 


Bo375 


SHP318.D1 


4.:! 


01) 


i() 


27.0 


12 


01 


51 


1490 


Bo:38() 


Slll'aiiJ 



" - source identifications beginning with Bo refer to Globular Cluster candidates listed in Table IV of Battistini et al. (1987), 

MIT - in Magnicr (1993). 

' - source identifications beginning with SHP refer to M31 ROSAT/PSPC X-ray source catalog entries from Supper et al. 
(1997, 2001). Identifications beginning with D refer to the list of GC X-ray sources listed in Table 3 of DiStefano et al. (2002). 
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Table 4 

Best-fit model parameters of the energy spectra of the bright GC sources. XMM-Newton/EPIC data, 0.3 — 10 
keV energy range. Absorbed simple power law model approximation. Parameter errors correspond to the la 

level. 



ID " Observation /Instrument 

Model: Absorbed Power Law (PO W ERLAVV ^ WAijIj) 





Photon 

Index 


< 
H 


Plux<= 


,X 
(d< 


i 

f) 


X 




2 


l,r,r, + 0,03 


0,1.' + 0,01 


28, «0 + 0,32 




LSI) 


1990 


xnnn 8 / { M f )H 1 -i- M f )S 2 ) 


(i 


1,^1 " ^' 

-0,10 


„■ 0,11,-. 


(KSl i,,Uli 






,-■0 




7 


1.70 + "--''-' 
-0.16 


U 31 + 0.00 
-0.06 


1.44 ± 0.07 


50.4 


49) 


99 


xmm 7/(MOSl+MOS2) 


8 


^•^^-0.15 

^■''''-0.18 
, ,„+0.08 
'^■^''-0.11 


o.45+0;0« 
„ 3.+0.02 

"■37_0.04 


0.93 ± 0.07 
0.80 ± 0.06 
1.65 ± 0.07 


30.2 
23.8 
72.6 


29) 
30) 
73) 


64 
55 
114 


xmm 1/pn 
xmm 4/pn 

xmm 6/pn 


9 


^•'*-^-0.28 


61 + "-^" 
-0.11 


0.58 ± 0.06 


18.2 


15) 


40 


xmm 1/pn 




2.44+0.20 
-0.26 
1 Rq+0.08 
^•^^-0.06 


0.54 ± 0.10 

o"tS:°o! 


0.74 ± 0.06 
3.25 ± 0.10 


29.7(33) 
83.5(105) 


51 
225 


xmm 4/pn 
xmm 6/pn 


10 


^•■^'^-0.07 


0.44 ± 0.03 


32.20 ± 0.64 


94.6 


94) 


2226 


xmm 4/pn 




1.52 ± 0.03 


o«lS:°o2 


31.40 ± 0.30 


231.2 


293) 


2170 


xmm 7/(MOSl+MOS2+pn) 


11 


1.42+""* 
-0.07 


0.09 ± 0.02 


9.72 ± 0.20 


115.2 


113) 


672 


xmm l/(MOSl+MOS2) 




1.34 ± 0.06 


0.09 ± 0.02 


9.45 ± 0.29 


72.5 


60) 


653 


xmm 2/pn 




1.61 ± 0.04 
1,56 ± 0,03 


o-14+^-^l 

— 0.04 
0.12 ± 0.01 


8.29 ± 0.25 
8.27 ± 0.09 


195.9 


189) 


573 
571 


xmm 4/(MOSl+pn) 
xmm 6/(MOSl + MOS2 + pn) 


14 


^ 90 + U.UV 
-0.06 


25+U.Ub 
-0.01 


1.87 ± 0.06 


125.7 


(97) 


129 


xmm l/(MOSl + MOS2 + pn) 




1.81 ± 0.06 


„.31 + 0.03 


2.20 ± 0.06 


157.7 


145) 


152 


xmm 4/(MOSl + MOS2 + pn) 




2.30+0.05 


031 + 0.04 


0.98 ± 0.04 


127.2 


138) 


68 


xmm 6/(MOSl + MOS2 + pn) 


16 


i-«*t^:o? 


o-i^:o.o3 


1.18 ± 0.05 


93.0 


89) 


82 


xmm l/(MOSl+MOS2+pn) 




i-^^+ois 


0.20+0-07 


2.73 ± 0.20 


34.6 


26) 


190 


xmm 2/pn 




1 91+00* 
— 05 


0.22+0 01 


3.36 ± 0.07 


184.5 


209) 


232 


xmm 4/(MOSl+MOS2+pn) 




2 17+0.04 
— 0.03 


0-27+0-01 
— 0.04 


2.12 ± 0.04 


208.8 


243) 


147 


xmm 6/(MOSl + MOS2 + pn) 


19 


l-^^+019 


< 0.06 


3.10 ± 0.28 


70.8 


62) 


21 


xmm 4/pn+xmm 6/pn 


20 


1.82+"-"7 
-0.08 


0.12 ± 0.02 


1 .99 ± 0.08 


238. 1 


214) 


137 


xmm 6 /pn 


22 


1.57+^-^^ 


0.26 ± 0.01 


52.95 ± 0.42 


443.7 


424) 


3651 


xmm 5/(MOSl + MOS2) 


24 


-0.04 


0.13 ± 0.01 


4.39 ± O.09 


213.4 


199) 


304 


xmm l/(MOSl+MOS2+pn) 




2.11 ± 0.11 


0.18+0-02 
— 0.05 


4.69 ± 0.19 


51.0 


63) 


324 


xmm 2/pn 




— 0.04 


0.14 ± 0.01 


5.21 ± 0.09 




284) 


360 


xmm 4/(MOSl+MOS2+pn) 




1 .92 i 0.02 


0.12 i 0.01 


5.68 i 0.08 


486.3 


450) 


393 


xmm 6/(]V[OS l + MOS2 + pn) 


25 


1.60 ± 0.04 


0.13 ± 0.01 


5.31 ± 0.10 


197.7 


201) 


367 


xmm l/(MOSl + MC:)S2 + pn) 




1.43+0.03 


0.09+0-02 


7.64 ± 0.11 


320.6 


316) 


528 


xmm 4/(MOSl + MOS2 + pn) 




1.57+0-02 
— 0.03 


0.13 ± 0.01 


7.15 ± 0.09 


524,7 


517) 


495 


xmm 6/(MOSl + MOS2 + pn) 


26 


o-t"o:"o'3 


0-o^troi 


8.91 ± 0.20 


228.0 


176) 


614 


xmm 5/(MOSl + MOS2) 
+xmm 6/pn 


28 


2.23+"/,^, 


0.12 ± 0.02 


1.94 ± 0.05 


81.7 


91) 


134 


xmm l/(MOSl + MOS2 + pn) 




2 25 + '' "-°^ 
-0,06 


0.12 ± 0.02 


1.87 ± 0.05 


132.7 


118) 


129 


xmm 4/(MOS2+pn) 




2.13 ± 0.06 


o-"iro? 


1.62 ± 0.04 


84.8 


87) 


112 


xmm 6/pn 


29 


2.12 ± 0.04 


0.14+0-0^ 
0.10+0:03 


2.80 ± 0.07 


145.4 


164) 


194 


xmm l/(MOSl+MOS2+pn) 




l.S2+°-°^ 
— 0.07 


4.52 ± 0.20 


36.0 


46) 


313 


xmm 2/pn 




1.61+003 


0.09 ± 0.01 


6.50 ± 0.11 


320.4 


333) 


450 


xmm 4/(MOSl+MOS2+pn) 




1.96 ± 0.04 


0.16+0-01 


4.19 ± 0.09 


315.0 


282) 


290 


xmm 6/(MOSl+MOS2+pn) 


30 


2 S7+0-^^ 
-0.28 


0-^troi 


0.35 ± 0.03 


34.9 


36) 


25 


xmm 4/pn 


31 


1.56 ± 6.02 


0.09 ± 6.01 


15.95 ± 0.18 


483.4 


463) 


1103 


xmm l/(M0gl+M0g2+pn) 




, ^„+0.04 
l''*-0.03 


0.10 ± 0.01 


16.56 ± 0.38 


163.0 


202) 


1144 


xmm 2/(MOSl+MOS2+pn) 




1 71 + 0.01 
-0,03 
1,66 ± 0,02 


0.10 ± 0.01 


13.31 ± 0.18 


554.7 


571) 


920 


xmm 4/(MOSl + MOS2 + pn) 




0.09 ± 0.01 


12.53 ± 0.16 


731.6 


759) 


866 


xmm 6/[MOSl + MOS2 + pn) 


32 


0.67 ± 0.04 

"--till 
0.56 ± 0.02 
0.57 ± 0.03 


0.12 ± 0.03 

"--till 

< 0.05 
0.04 ± 0.02 


18.64 ± 0.38 

20.03 ± 0.66 

26.07 ± 0.34 
10.82 ± 0.19 


207.0 
90.6 
412.4 
326.0 


182) 

48) 

348) 
249) 


1288 
1384 
1802 
748 


xmm l/(MOSl + MOS2 + pn) 

xmm 2/pn 
xmm 4/(MOSl+MOS2+pn) 
xmm 6/(MOSl+MOS2+pn) 


34 


-«t"o:f3 


o-otro'4 


0.27 ± 0.02 


53.1 


62) 


19 


xmm 1/pn+xmm 6/pn 


36 


-al^o:f6 


0.23+0-^ 


0.66 ± 0.07 


20.7 


33) 


46 


xmm 6/pn 


37 


--i^o-.l!^ 


0.07 ± 0.01 


10.87 ± 0.17 


181.8 


200) 


751 


xmm l/(MOSl+MOS2+pn) 




1.55+0-10 
-0,09 


0.08 ± 0.02 


8.77 ± 0.34 


33.8 


39) 


606 


xmm 2/pn 




1 44 + 0-03 
-0.02 


0.09 ± 0.01 


10.85 ± 0.15 


317.1 


282) 


750 


xmm 4/(MOSl + MOS2+pn) 




1.68 ± 0.03 


0.12 + 001 
-0.02 


9.24 ± 0.10 


401.7 


352) 


640 


xmm 6/(MOSl + MOS2+pn) 


38 




0.23 ± 0.05 


0.75 ± 0.05 


33.2 


32) 


52 


xmm 5/pn 


39 




0.04 ± 0.03 


0.64 ± 0.03 


26.3 


26) 


44 


xmm 5/(MOSl + MOS2) 


40 


1.90 ± 0.32 


0.19 ± 0.07 


0.54 ± 0.05 


14.6 


16) 


37 


xmm 6/pn 


43 


1.56 ± 0.04 


0.13 ± 0.01 


21.80 ± 0.30 


238.3 


199) 


1507 


xmm 9/(MOSl+MOS2+pn) 



" - Source number in Tabic 3 

^ - Equivalent hydrogen column depth in units of 10^^ cm"^ 

— Absorbed model flux in the 0.3 — 10 keV energy range in units of lO"-*^^ erg s"-*^ cm"^ 

— Absorbed isotropic source luminosity in the 0.3 — 10.0 keV energy range in units of 10^^ erg s"-*^ assuming the distance of 760 kpc 



15 



Table 5 

Best-fit model parameters of the energy spectra of the bright GC sources. Chandra/ ACIS data, 0.5 — 7.0 keV 

ENERGY RANGE. ABSORBED SIMPLE POWER LAW MODEL APPROXIMATION. PARAMETER ERRORS CORRESPOND TO THE la LEVEL. 



ID " Observation/Instrument 
Model: Absorbed Power Law (POWERLAW*WABS) 
Photon Fh^ P If 



Index (dof) 



G 


2.0 il;;:- 


< 0.08 


0.19^0.01 




17) 


:il 


2()i7/AC'lS 


7 


±.»0_Q 25 


U.DU_Q 20 


0.73 ±0.06 


31.2 


.23) 


51 


2017/ACIS 


8 






'J-'J^J-o.os 


1.15 ±0.06 


21.3 


;26) 


80 


1575/ACIS 


9 


-0.09 


'^•^i-o.oe 


2.07 ±0.08 


20.9 


;39) 


143 


1575/ACIS 


10 


1.47 ±0.15 


0.56 ±0.14 


25.35 ± 1.00 


27.4 


.28) 


1752 


1580/ACIS 




1.481°:;! 


0.641?-° 


28.48 ± 1.07 


27.6 


;25) 


1968 


1584/ACIS 




1-461H! 


0.49 ± 0.03 


33.72 ± 0.40 


238.6 


;i76) 


2330 


2017/ACIS 




1 49+011 

1^-^^-0.13 


'-'■'J^-0.11 


28.57 ± 1.05 


31.8 


;32) 


1975 


2901/ACIS 


11 


1.44 ± 0.06 


0.14 ±0.03 


7.87 ±0.18 


86.4 


;65) 


544 


1575/ACIS 


14 




Q 01 +().uy 

'^•'Jl-0.()4 


0.91 ± 0.05 


19.6 


;i5) 


63 


1575/ACIS 


16 


1- ' '^-0.09 


24+"''" 


3.94 ±0.14 


37.2 


;34) 


272 


1575/ACIS 


20 


1.9011!:^^ 


o.2o;;;:i| 


0.50 ±0.04 


10.8 


;io) 


35 


1575/ACIS 


22 


1.48 ± 0.08 
1.661°-;? 
1.51 ±0.03 


n 99+U.U3 
^ — 02 

0.35l°-?« 
0.30 ± 0.02 


49.48 ± 1.23 
47.61 ± 1.46 

59.80 ±0.56 


60.8 
34.8 

261.7 


;62) 

'38) 

'215) 


3420 
3290 

4130 


309/ACIS 
1582/ACIS 

1575/ACIS 


23 


1- ' 1-0.12 


< 0.06 


0.54 ±0.04 


14.6 


'19) 


38 


1575/ACIS 


24 


1.84 ±0.05 


1 5+°-°^ 

U-lO-O.Ol 


6.02 ±0.12 


90.3 


;84) 


416 


1575/ACIS 


25 


1.36 ± 0.05 


0.07 ±0.01 


7.87 ±0.15 


78.8 


;84) 


544 


1575/ACIS 


26 


9l+° "7 


'J-^^-o.os 


8.37 ±0.25 


76.4 


;6i) 


578 


1575/ACIS 


27 


2.02in^ 


0.04 ± 0.02 


1.42 ±0.05 


33.0 


;25) 


98 


1575/ACIS 


29 


1.841-- 
1.96 ± 0.07 


^•l'5-O.Ol 


2.64 ±0.16 

3.65 ± 0.09 


20.4 
76.6 


;25) 
;56) 


183 
252 


303/ ACIS 
1575/ACIS 


30 


1 69+"- ^"^ 

±.U»_(, 27 


'^•^"-0.09 


0.43 ± 0.04 


12.6 


;i6) 


30 


1575/ACIS 


31 


1--*1-0.04 

1.52 ±0.04 


0.11 ±0.03 
0.06 ±0.01 


15.88 ± 0.44 
14.21 ±0.20 


54.7 
210.4 


;59) 

'171) 


1097 
982 


303/ACIS 
1575/ACIS 


32 


o.eolil:;^ 

0.26 ±0.14 


< 0.08 
03+° " 

'-'•^"^-0.03 


16.27 ±0.99 
27.23 ± 1.52 


18.3 
10.3 


;2o) 
;i9) 


1125 
1880 


310/ACIS 
1854/ACIS 




o.5oi;::;t 


0.261°:°^ 

o.isl":- 


8.68 ±0.34 
14.12 ±0.67 


82.1 
47.2 


'85) 
;33) 


600 
976 


2017/ACIS 
2494/ ACIS 


35 




0.11 ±0.04 


15.30 ±0.71 


30.9 


;23) 


1057 


1582/ACIS 


37 


1.40 ± 0.04 


< 0.04 


10.32 ± 0.23 


102.8 


(93) 


713 


1575/ACIS 


41 


1 44+U-i'> 

1 cc+o.ig 

l.dd_(, 14 


< 0.09 

0.041°:°? 


16.02 ± 0.69 

16.84 ±0.88 


29.3 

10.9 


'28) 

;i7) 


1107 

1164 


2895/ACIS 

2896/ACIS 






0.081°;°^ 


16.35 ±0.77 


22.2 


;22) 


1130 


2897/ACIS 


42 


1.66 ±0.06 


0.28 ±0.03 


74.50 ± 1.43 


134.2 


;io8) 


5148 


2052/ ACIS 




1 yr + O.OT 
l-'^-0.06 


n 90+0.03 

0-2»_o.02 


150.0 ± 3.0 


102.9(69) 


10372 


2053/ACIS 



" - Source number in Table 3 

' - Equivalent hydrogen column depth in units of 10^^ cm~^ 

° - Absorbed model flux in the 0.3 — 10 keV energy range in units of 10~^^ erg s~^ cm~^ 

- Absorbed isotropic source luminosity in the 0.3 — 10.0 keV energy range in units of 10^^ erg s~^ assuming the distance of 
760 kpc 
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Table 6 

Bright GCS spectral fit results {XMM-Newton/EPIC data, 0.3 — 10 keV energy range; Chandra/ ACIS data, 
0.5 — 7.0 keV energy range). Absorbed cut-off power law and Comptonization models. Parameter errors 

CORRESPOND TO la LEVEL. 



ID " 
















Remarks 


Model: Absorbed Cutoff Power Law (CUTOFFPL*WABS) 




Photon 


Cutoff Energy 






Flux= 




Td 






Index 


(keV) 








(dof) 






22 


o.64t^:i^ 






16+""^ 


53.07 ±0.52 


228.2(215) 


3668 


Chandra 1575 










16+0"2 


47.74 ± 0.40 


399.3(423) 


3300 


xmm. 5 


26 


0.18 ±0.09 


4.861^-^2 




< 0.02 


7.91 ±0.17 


208.0(175) 


545 


xrnm 5+xmni 6 


32 


-0.15 ±0.09 






< 0.03 


16.71 ± 0.28 


171.8(181) 


1155 


xmm 1 




0.17 ±0.15 


K 4Q+2.05 




< 0.02 


18.59 ± 0.61 


77.4(47) 


1285 


xmm 2 




0.21 ±0.07 


7 +1.76 




< 0.01 


23.71 ± 0.31 


376.6(347) 


1640 


xmm 4 






r 74+0.66 
^•'^-1.05 




< 0.01 


10.69 ±0.17 


278.0(248) 


739 


xmm 6 


42 




■J-OJ_0.97 




0.16 ±0.05 


68.39 ± 1.32 


126.4(107) 


4726 


Chandra 2052 




1 3+0-31 


l-8li°o1f 




0.09 ± 0.03 


125.9 ± 2.6 


74.2(68) 


8700 


Chandra 2053 


43 


1 13+"-^''' 


6.82tt:SI 




0.09 ± 0.02 


20.39 ± 0.29 


232.0(198) 


1409 


xmm 9 


Model: Absorbed Comptonization Model (COMPTT*WABS) 




fcT„ 




r 




Flux= 


x' 








(kcV) 


(keV) 








(dof) 






22 


20+" "^ 


J--'*O_0.07 


26.3+2^ 


<J-J-'_0.07 


50.29 ± 0.49 


229.2(217) 


3475 


Chandra 1575 




'-'■'J' -0.03 


1 OO+0.23 
i.O»_o.i8 


19.5i;j 


< 0.09 


47.82 ± 0.41 


398.7(422) 


3297 


xmm 5 


26 


03+" "'^ 


^- '"-0.12 


34.6+^2-5 




7.38 ± 0.16 


206.8(174) 


509 


xmm 5+xmm 6 


32 


0-06Io.o6 
07+° °'^ 
09+° °^ 

in+° °'5 

'^•J^'J-o.og 


J^-o'-o.io 
1 63+013 

1.00_Q jg 

1 09+0. 08 
■L.O^_0.06 

l- ' ="-0.07 


38.6lf, 
44.8ti» 
43.211;? 
45.21^;^ 


< 0.09 

< 0.06 


16.11 ± 0.27 
17.63 ± 0.58 
21.63 ± 0.28 

9.93 ±0.16 


177.0(180) 
65.5(46) 
341.7(346) 

252.3(247) 


1113 
1220 
1495 

686 


xmm 1 
xmm 2 
xmm 4 

xmm 6 


42 


"•J^'-0.()6 




23.7+f^ 


19+0-03 

'-'•-'-'^-0.06 


65.46 ± 1.26 


125.7(106) 


4524 


Chandra 2052 




< 0.14 




29.51;:;;; 


0.18 ±0.03 


122.8 ± 2.5 


71.2(67) 


8485 


Chandra 2053 


43 


0.14 ±0.04 




22.4t^;^ 


0.09 ± 0.02 


19.42 ± 0.28 


258.3(235) 


1339 


xmm 9 



" - Source number in Table 3 

' - Equivalent hydrogen column depth in units of 10^^ cm~^ 

° - Absorbed model flux in the 0.3 — 10 keV energy range in units of 10"^^ erg s"^ cm~^ 

- Absorbed isotropic source luminosity in the 0.3 — 10.0 keV energy range in units of 10*® erg s"^ assuming the distance of 
760 kpc 
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Table 7 

Bright GCS spectral fit results, XMM-Newton/EPlC and Chandra/ ACIS data. Two-component model 
approximation: (BB0DYRAD+DISKBB)*WABS and (BB0DYRAD+P0WERLAW)*WABS models. Parameter 

errors correspond to la level. 



ID " 
















Observation 


Model: (BBODYRAD+DISKBB)*WABS 




kTBB 

(keV) 


kTin 

(keV) 


rinVcosi 
(km) 




FIux= 


(dof) 


rd 
J-'X 




22 


1.40i^;^^ 


0.83 ±0.16 


37t^^ 


0.17 ±0.03 


50.94 ± 0.50 


229.3(217) 


3520 


Chandra 1575 






J^-"0-0.14 


26l| 
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Fig. 1. — Optical image of M31 from the Digitized Sky Survey with regions covered by XMM-Newton/EPlC {black circles) and Chan- 
dra/ ACIS/BRC observations (white regions). The remote globular cluster Gl field is located outside image boundary. 
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Fig. 2. — Combined XMM/EPIC images of M31 fields in the 0.3 — 7.0 keV energy band. The image was convolved with circularly symmetric 
Gaussian function with a = 8" . The globular cluster X-ray sources detected in XMM-Newton and Chandra observations are marked with 
circles. 
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Fig. 3. — Combined Chandra/ACIS images of M31 fields in the 0.3 — 7.0 kcV energy band. The image was convolved with circularly 
symmetric Gaussian function with u = 8". The globular cluster X-ray sources detected with Chandra are marked with circles. 
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Fig. 4. — Upper panels: Representative XMM-Newton/EPIC spectra of the bright GC X-ray sources in the 0.3 — 10 kcV energy band. 
(EPIC-pn is always the upper spectrum and the two MOS spectra overlap). Corresponding model fits to the EPIC-pn, MOSl and MOS2 
spectra (Tables 4, 6, 7) are shown with thick, thin and dotted histograms respectively. Left panel: Bo 96 (#14), EPIC-pn and MOSl, 2 
data, fit by absorbed simple power law model. Middle panel: Bo 135 (#10), EPIC-MOSl and MOS2 data, fit by absorbed power law 
model with exponential cutoff. Right panel: Bo 158 (#32), EPIC-pn and MOSl, 2 data, fit by absorbed Comptonization model. Lower 
panels: Chandm/ ACIS spectra of the bright GC X-ray sources. Model fits arc shown with thick histograms (Tables 5, 6, 7). Left panel: 
Bo 78 (#9), ACIS-S data, 0.5 — 7 kcV energy range, fit by an absorbed simple power law model. Middle panel: Bo 82 (#10), ACIS-S 
data, 0.5 — 10 keV energy range, fit by an absorbed simple power law model. Right panel: Bo 375 (#42), ACIS-S data, 0.5 — 10 keV 
energy range, fit by an absorbed power law model with exponential cutoff. 
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Fig. 5. — The distribution of the spectral slopes derived from the analysis of 31 bright GC X-ray sources. The histogram bins are 0.1 wide. 



23 



I I I I I 1 1 — I — I I I I I I 1 1 — I — I — I I I I I 1 1 — I — I — I I I I 



2.5 - 



X 
I — f 

c 
o 



o 

I 1 

-I-' 

o 

CD 



1.5 



1 - 



M31 Globular Cluster Sources 



-f 



0.5 



I I I I I I I I I I I I I I I I I I I I I I I I I I 



10 100 1000 

X-ray Luminosity, xlO^s ergs/s 



10000 



Fig. 6. — Spectral photon index of the GC X-ray sources in M31 vs. their X-ray luminosity in the 0.3 — 10 keV energy band. The error bars 
in X-axis reflect statistical uncertainty of the source flux determination and in some cases the range of source X-ray luminosities observed 
with XMM and Chandra. 
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Fig. 7. — Left panel: The values of A'^H derived from X-ray spectral fitting plotted against the values of A^h derived from optical extinction 
data. The dotted line shows points with equal X-ray and optical A^h- The dashed lines mark expected Galactic foreground absorbing column 
in the direction of M31. Right panel: The distribution of the absorbing columns derived from the spectral analysis of 31 bright GC X-ray 
sources. Each bin along X-axis has a width of 5 X 10^" cm~^. The Galactic foreground absorbing column in the direction of M31 (7 X 10^'' 
cm~^) is marked with dotted line. 
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Fig. 8. — X-ray flux histories of 10 GC sources in our sample obtained combining tlic data of Chandra/ ACIS {filled circles), HRC {open 
circles) and XMM/EPIC {open rectangles) observations, 0.3 — 10 keV energy range. Tlic upper limits correspond to a 2a level. 
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Fig. 9. — X-ray flux histories of the recurrent M31 GC source Bo 163 (#35) based on XMM-Newton and Chandra {upper panel) and ROSAT 
{lower panel) archival data. The upper limits correspond to a 2cr level. 
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Fig. 10. — Spectral variability of M31 GC sources Bo 78 (#9)(ieft panel) and Bo 148 (#29) {right panel). The hardness of the spectrum 
expressed in terms of the spectral photon index (Y-axis) is plotted against source X-ray luminosity in the 0.3 — 10 keV energy band (X-axis). 
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Fig. 11.— Upper panel: X-ray light curve of Bo 86 (#11) during the 2002 January 6 XMM-Newton/EPIC-pn observation (0.3 - 10 keV 
energy range and a 1000 s time resolution). Middle panel: X-ray light curve of Bo 107 (#16) during the 2000 June 25 XMM-Newton/EPlC- 
pn (0.3 — 10 IjeV energy range and a 1000 s time resolution). Lower panel: X-ray light curve of Bo 135 (#22) during the 2001 October 5 
Chandra/ ACIS observation (0.3 — 7 keV energy range and a 2000 s time resolution). 
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Fig. 12. Broad-band color (hardness) vs. intensity for three XMM-Newton observations of X-ray souree in the globular cluster Bo 148 
(#29). The broad- band color is defined as ratio of source intensities in the 2.0 — 7.0 and 2 — 7 keV energy bands. The EPIC-MOS data is 
binned to 1500 s. The corresponding source luminosity changes between ~ lO^'^ and ~ 4 x lO^'^ ergs s~^ in the 0.3 — 10 keV energy band. 
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Fig. 13. — Upper panel: X-ray light curve of Bo 158 (#32) during the 2002 January 6 XMM-Newton observation, obtained from combined 
EPIC-pn, MOSl, and MOS2 cameras (0.3 — 10 keV energy range and a 500 s time resolution). Middle panel: X-ray light curve of Bo 158 
during the 2001 July 24 Chandra/ ACIS observation (0.3 — 7 keV energy range and a 1000 s time resolution). Lower panel: The X-ray 
modulation fraction (%) plotted as a function of X-ray luminosity of Bo 158 (in units of 10^* ergs s~^) 
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Fig. 14. — Cumulative X-ray luminosity distributions of M31 GC X-ray sources and central bulge of M31. 
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Fig. 15. — Spatial distribution of M31 GC X-ray sources detected with XMM-Newton and Chandra. The brightest sources with luminosities 
above 10*^ ergs are marked with large white circles. X-ray sources with luminosities occasionally exceeding 10^® ergs are shown with 
white boxes. The positions of fainter GC X-ray sources are shown with small white circles. 
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Fig. 16. - - Left panel: V-inagnitudc distributions for M31 GC candidates. The results for GC hosting bright X-ray sources are shown with 
thick histogram. The distribution for the GC candidate sample from Battistini et al. 1987 is shown with dotted histogram. Right panel: 
metallicity ([Fe/H]) distribution for M31 GC candidates hosting bright X-ray sources. 
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Fig. 17. — Left panel: the measured isotropic 0.3 — 10 keV X-ray luminosity of M31 GC X-ray sources versus the host globular cluster 
metallicity ([Fe/H]). The error bars in Y-axis reflect statistical uncertainty of the source flux determination and in some cases the range of 
source X-ray luminosities observed with XMM and Chandra. Right panel: the spectral power law photon index for the GC sources from our 
sample vs. metallicity of the globular clusters hosting them. 



